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Abstract. We study the solutions ol Von Neumann's expanding model with 
reversible processes for an infinite reaction network. We show that, contrary to 
the irreversible case, the solution space need not be convex in contracting phases 
(i.e. phases where the concentrations of reagents necessarily decrease over time). At 
optimality, this implies that, while multiple dynamical paths of global contraction 
exist, optimal expansion is achieved by a unique time evolution of reaction fluxes. This 
scenario is investigated in a statistical mechanics framework by a replica symmetric 
theory. The transition from a non-convex to a convex solution space, which turns out 
to be well described by a phenomenological order parameter (the fraction of unused 
reversible reactions) is analyzed numerically. 
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1. Introduction 

At the most basic level of abstraction, a transformation process like an industrial 
production technology or a chemical reaction can be specified merely by the coefficients 
measuring the amounts of each commodity or chemical species that are consumed and, 
respectively, produced when the process operates at, say, unit scale (or flux). When N 
such processes are brought together in a network connecting M species so that the inputs 
of one process may be the outputs of another, the emergent properties of the network will 
depend crucially on the particular choice of the (quenched) "stoichiometric" coefficients. 
Of special interest in this respect are questions regarding network optimality. In a 
landmark paper [I], J. Von Neumann considered the problem of finding the largest 
(uniform) species production rate possible for irreversible processes specified by given 
stoichiometry. The existence of a maximal expansion rate p* can be proven easily 
for any matrices of input-output stoichiometric coefficients satisfying broad generic 
assumptions [2]. 

The study of this problem with random stoichiometry [3j[4] provides hints on the 
typical behavior of the solution of the Von Neumann problem in complex situations. In 
the limit where the number (N) of reactions and (M) of species diverge, with a fixed 
ratio {n = N/M) a full characterization of the growth properties is possible, using tools 
of statistical physics. The maximal growth rate is given, to a first approximation, by 
the ratio of output and input coefficients, with a non-trivial correction which is of the 
order of l/\/K, where K is the (average) number of reactions in which each species 
is involved. This latter term depends on the structure of the network of reactions. In 
particular, as the ratio of processes to species increases, the maximal growth rate also 
increases. In addition, one can show that the space of solutions for a given growth rate 
is convex and that it shrinks to a single point when the growth rate approaches the 
maximal one, for a specific network of reactions [3j[4]. 

When the network of reactions is constrained to obey the mass conservation implicit 
in chemical reactions [H], expanding solutions are not possible, as the maximal growth 
rate cannot be larger than zero. The Von Neumann problem then reduces to that of 
metabolite producibility in metabolic networks, addressed in |6|7j. The analysis of fluxes 
in the metabolic network of the bacterium E. coli within Von Neumann's framework 
has been discussed in [8] . Interestingly, the presence of conserved pools of metabolites 
imply that the solution space for the maximal growth rate does not coincide with a 
single point. Rather, a finite volume of flux configurations corresponding to maximal 
growth rate (equal to zero) exists. We refer the interested reader to [8) for the biological 
implications. Here we merely point out that even in this case the space of solutions 
retains its convexity properties: any linear combination of two solutions with positive 
weights is still a solution. 

These results apply to a network of irreversible reactions (or generically when 
reversible processes can be split into two separate reactions). This paper addresses the 
case where one or more of the reactions can also be run in the reverse direction. Such a 
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scenario is not usually considered in economic applications, as technologies are generally 
unidirectional. For biochemical systems it is however crucial, both because many of the 
chemical reactions occurring in cells are physiologically reversible and because the shape 
of the solution space, particularly its convexity properties, has strong implications on 
the effectiveness of algorithms designed to find the optimal flux vectors and sample the 
solution space uniformly. 

We first show how the case in which processes can be reversed is a non-trivial 
generalization of the irreversible case. In particular, convexity is no longer guaranteed 
when the growth rate is negative, whereas it can still be proven when the growth rate 
is non-negative. In order to address what happens in typical cases, we replicate the 
analysis on random systems. Reversibility turns out to generate a substantially more 
complex scenario for the solution space of Von Neumann's expansion problem, which 
considerably complicates its numerical analysis. For simplicity, we focus on the fully- 
connected version, defined in Sec. 2, resorting to the replica trick to compute the 
maximal growth rate (Sec. 3). In Section 4 the analytic solution is discussed and 
compared with numerical results obtained for a large network with a given stoichiometry. 
This analysis confirms that, when the growth rate is negative, the space of solutions 
splits into many disjoint components. A brief outlook is given in the concluding Section 
5. Finally, an Appendix details the heuristic algorithm employed for the numerical 
exploration of the solution space. 

2. The problem and a basic observation 

Let us consider N processes, labeled by i = 1, . . . , N, which operate on a set of 
M species - be them chemical substances or commodities - labeled by p = 1, . . . , M. 
These are transformation processes whereby some inputs are transformed into some 
output products in proportions which are specified by the stoichiometric matrices of 
outputs A = {af} and inputs B = {fof}. Time is discrete, and we consider the case 
where all processes are run in parallel at scale o~i(t) > 0, between time t and time t+1. 
This means that <Tj(i) processes of type i are run in parallel, consuming o~i(t)b^ units of 
species p = 1, . . . , M and producing £Xj(t)af units of species p' = 1, . . . , M. 

Hence, the total quantity of species p consumed as input at time t, by all processes, 
is = J2iLi a i(t)K an d the total amount of output of species p produced is 

M (t) = Sili a i{t) a i ■ Given an initial concentration O M (0) of species, feasible growth 
paths {cr(t)} t> o (with cr{t) = {o" i (t)}^ 1 ) are those such that cjj(t) > for all i and t, and 
that the process at time t produces enough outputs to run the processes at time t + 1, 
i.e. 0^(t) > I^{t + 1) for all t. Von Neumann [l] further focuses on paths o~i{t) = Sjp*, 
with Si > 0, of exponential expansion (p > 1) or contraction (p < 1). Feasible solutions 
s with growth factor p (and growth rate log p) are easily seen to satisfy the constraints 

N 

c"(s) = s i « - PK) > 0, Vp = 1, . . . , M (1) 

i=l 
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or, more compactly, (A — pB)s > 0. If p is very small, the space of solutions is very 
large, in particular any s > is a solution for p = 0. On the contrary, if p is large 
enough no vector s will satisfy the conditions 0. For any A and B satisfying broad 
generic assumptions [2] it is easy to show that there is a maximal expansion rate p*. 
Von Neumann's (VN) problem amounts to finding the maximal growth rate and the 
path of maximal growth, i.e. 

max p subject to (A — pB) s > 0. (2) 

s>0 

Let us now generalize this problem to introduce reversible processes. Reversibility 
implies that input and output coefficients can be interchanged. If process i is reversible, 
it is possible to produce a quantity bf of species p by consuming units of species 
p! . It is interesting to observe at the outset that reversing the direction of reactions 
does not correspond to merely reversing the direction of time. Indeed, imagine that a 
solution s with growth rate p exists. Then one would naively expect that, if all reactions 
are reversed, the situation where all processes are run in the same proportions (i.e. with 
operation scales s = s) should yield a solution with p = 1/p. Indeed, it is not difficult 
to find that the maximal growth rate p corresponding to the reversed solution satisfies 

^HK 1 -^) • o *< s >=f>- < 3) 

which means that in order to obtain a growth rate p = 1/p for the reversed process the 
condition d 1 = for all p is necessary but not sufficient]!} Hence reversing reactions is 
not, in general, equivalent to time reversal. 

By the same arguments as above, one finds that, when reactions can be run in both 
directions, Von Neumann's (VN) problem can be cast as 

maxp subject to S p (s)s > (4) 

s 

where s is a flux vector and H p (s) is a M x N matrix with entries given by 

{<-"»? fOTS .>° (5) 
I bf — pa% for Si < 

Clearly, because the matrices depends on scales s, the above conditions are not linear in 
s. This generates key differences between the irreversible case and the present model. 
In the former case, in fact, a convex combination As + A's' of two solutions s and s' with 
non negative coefficients A and A' is always a feasible solution of the original problem. 
Irreversibility can modify this picture. To see this it is convenient to introduce the 
quantities 



N 



c » = V - pK) - e(-sm - p<)\ (6) 



i=l 



| In order to show this, let c M (s) = J^i s i{^i ~ P a t) an d note that 

c^(s)+ pc^(s) = (1 - /»p)0"(s) 

This immediately tells us that if c M (s) = c M (s) = for all /i, then p = 1/p. If however c M (s) > 0, one 
has (^(s) < (1 — pp)O ll {s). Eq. (|3j) is derived straightforwardly from this. 
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which are required to be non-negative for each p under conditions Q. Indeed 
substituting one finds 

c"(As + A's') = Ac^(s) + A'c^(s') + (1 - p)^(s, s') (7) 

where 

•A"(s, s') = J2« + K)[^9(~s t ) + AV^C-sJ) - (As, + AVJ^-A* - A's-)] (8) 

i 

If s and s' are both solutions, then a sufficient condition for their linear combination to be 
again a solution is that (1 — p)A^(s, s') > 0. If Sj and have the same sign, for all i (i.e. 
Sis'i > 0), then ^4 M (s,s') = 0, which means that convexity is guaranteed for all values of 
p. But if the two solutions differ by a change of sign (i.e. if SjS^ < for some i), then the 
term in square brackets in Eq. (j8l) is non positive and it equals — min[A|sj|, A'|s^|] < 0. 
This implies that convexity (i.e. (1 — p)A tl (s, s') > 0) is guaranteed only if p > 1. Note 
that convexity may not occur even at p* when p* < 1. 

Hence when p < 1 it is possible that the space of feasible solutions is disjoint 
in different connected pieces. However the condition (1 — p)A fl (s,s') > for all p 
is sufficient but not necessary, as c^(As + A's') > may be satisfied for all p, even 
if the former condition is not satisfied. A trivial case when this happens is the case 
p = 0, where clearly all s e IR are solutions. Hence, in order to understand whether 
non-convexity is realized in typical cases, we're going to study ensembles of random 
instances in the next section. 



3. Typical properties of feasible solutions: Replica-symmetric theory 

In what follows we are going to take 

< = *(l + 4) (9) 



K = 5 (i + §) do) 

where af and /3f are going to be drawn, independently from some distribution with finite 
first and second moments. We assume that <pN reactions are reversible and (1 — <j))N 
are not (0 < < 1). To leading order in N we have p* = max{|, |} if > 0, which can 
be achieved by running the reversible reactions in the maximal growth direction while 
switching off (sj = 0) all reactions with growth rate min{|, |}. The nontrivial part of 
the growth rate is hence related to terms of order l/\fN and, in order to focus on those, 
we shall set a = b. As in the irreversible case [3], to account for subleading effects we 
set 

^ 1 + 7N 



and focus on g. The maximum allowed g shall be denoted by g*. 
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When N, M — > oo with n = N/M fixed, one can apply the replica trick to calculate 
g* (which is expected to be a self- averaging quantity). We can argue as in [3j: the 
volume of solutions for a specific choice of the stoichiometry is 

M / N \ 

V(p)=Tr s He[c»( S )]8[J2\si\-N) (12) 

M=l \i=l / 

where the trace over reversible reactions involves integrals from — oo to +oo, that over 
irreversible ones from to +oo. The 5-function enforces the constraint Yli l s «l = ^> 
which sets a scale for the fluxes. In the thermodynamic limit one expects the typical 
volume of solutions to be given by V oc exp(iVt>(p)) where 



v(p) = lim —\ogV(p) , (13) 

A 7 — >oo iV 

the over-bar denoting an average over the quenched disorder By the replica 

trick 



v{p) = lim lim — log V(pY (14) 

7V-s>oo r->0 J\ r 

After expressing the 9 function via its Fourier decomposition, carrying out the disorder 
average and isolating the emergent order parameter 

1 N 



N 



one arrives at 



V(pY = / Ji(q) J 2 (q)dq (16) 



/oo poo j£ 

Dz / DcJJexp 
Jo j 



AT \ /AT 



(17) 



J 2 (q) = Tr s 5 ^ |s« | - ^ Ml 5 S ~ ^ 



,i=l / \i=l 



where £, £' run from 1 to r and k = (af — /3f ) 2 . 

Under the replica symmetric (RS) Ansatz where 

<?«' = q + X&W (19) 
one finds that at the relevant saddle point 

v(g) = max^! + maxF 2 ] (20) 

<1,X rn,P,T 
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where m, (3 and r are additional order parameters, derived from imposing RS Ansatze 



on the Lagrange multipliers enforcing the 5-functions in (18), and 

(c + g + i^kq) 2 



H x = -E { log 



dc 



cxp 



2k X 







+ </>E € log 



2 

ds exp 



X + (1 — 0)E^ log / cis exp 



/5 2 / 

--s -(m 



/5 2 || 

s — m\s\ — rsf 

2 



where 



-« 2 /2 



(21) 

(22) 
(23) 



First note that for k = the problem has the trivial solution g = 0. Indeed, in the 
limit A; — ?- the integral in H\ is dominated by c = 0, and its value is largest for g = 0. 
So it is precisely the presence of fluctuations in input-output coefficients that allows for 
expanding phases or imposes contracting phases. 

The last two terms in H2 contain contributions of reversible and irreversible 
processes, respectively. The last term arises from the integration on the variable G M 
of a representative reversible process, whereas in the last but first term the integration 



is limited on Sj > 0, as appropriate for reversible processes. The problem in Eq. (20) 
has to be solved numerically for each values of the parameters. 



3.2. Optimal growth solution 

The solution simplifies in the limit p — > p* (or g — > g*) of maximal growth rate. 
Indeed, assuming that a unique flux vector survives in this limit, amounts to studying 
the solutions in the limit \ — > 0, since \ is proportional to the Euclidean distance 
between two solutions and S£i: 



(24) 



c*(0 



Evaluating H\ in this limit one finds that the integral is dominated by values of the 
constraints c given by 

otherwise 
so that for the distribution of c one finds 

p{c) 



2qk 



-9(c) + c 5(c) 



11.3* 
— — erf . : . 
2 2 ^Wq 



(25) 



(26) 



^/l-aqk 

For fluxes, it is convenient to study reversible and irreversible processes separately. In 
order to extrapolate the dominant contribution in the limit x ~ >* 0, we set 

b = fix , t = T X , tz= -mx- (27) 
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In the case of irreversible fluxes, one has 

f t(z - 0/6 for f < z 

C(0 = 1 (28) 

I otherwise 
and the corresponding distribution of fluxes is given by 

_ (s-tz/b) 2 

0e 2(t/t>) 2 ]_ ^ 

Pirr(s) = ^0,irr5(s) + 9(s) — , ^ ,irr = ~erfc^= (29) 

Some greater care is required for reversible reactions. It turns out that, at the optimal 
growth p*, the integral in s for reversible fluxes is dominated by 

t(z — £)/b for £ < min{0, z} 

s*rev(0 = { t(-z - 0/b for £ > max{0, -z) . (30) 
otherwise 

Note that for z > the support of the distribution of fluxes excludes the interval 
[—tz/b,tz/b] around the origin. Therefore, in such conditions solutions with null 
reversible fluxes are forbidden, i.e. each reversible reaction is active either in one 
direction or the other. If z < 0, instead, the support of the distribution extends over 
the whole real axis, but a finite fraction ipo >r ev of reversible fluxes are zero. Finally, one 
finds 



where 



Prev{s) = 1p ,revS(s) + 9{z)f + {z) + 6(-z)f.{z) (31) 



b e 2(t/f>) 2 be 2 <*/ fc ) 2 

'♦w-'M-^+^-ivir- (32) 

(s-tz/b) 2 (s + tz/b) 2 

be 2(t/6)2 be 2 ('/f») 2 

/_ (z) = 6l( a + tz/b) +e(-s + tz/b) (33) 



while 

V'o.ret; — Probj^ G [min{0, z}, max{0, — z}] | (34) 

is the probability that a reversible process is inactive (to which only negative values of 
z contribute). Summing up, the distribution of fluxes is given by: 

p( S ) = (pp rev (s) + (1 - (p)Pirr(s) (35) 

All of the above quantities can be evaluated directly once saddle point equations are 
solved for the order parameters and for g* as a function of n and 0. 



4. Results 



The numerical solution of the saddle point equations, in the optimal growth 
case, exhibits the following scenario. An expanding phase with g* > is possible if 
n > 1/(1 + 0) or n e f = n(l + 0) > 1, while the system is necessarily confined to a 
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0.1 1 10 

n(l-Kj)) 

Figure 1. g* j\fk versus n e / = n(l + </>). Lines represent analytic (replica) predictions, 
markers (with bars corresponding to the numerical error) results from simulations. 
Inset: detail for n e j < 1. 

contracting regime if n e j < 1 (see Fig. Q. This is in line with the results derived for 
the irreversible case, in view of the fact that N(l + <f>) is the total number of processes 
available to the system. Increasing 0, at constant n e f, i.e. increasing the share of 
reversible processes, appears to have a beneficial effect in the contracting regime (larger 
(J) means larger g* ), while the opposite occurs in the expanding regime (see Fig. [T]). 

At g*, the fraction of reagents that are not produced ('intermediates', with & 1 = 0) 
behaves similarly to the purely irreversible model as a function of n e f (see Fig. |2j left). 
In brief, as higher growth rates become achievable, the process becomes more efficient, 
with a decreasing fraction of "wasted" reagents (c M > 0). For the fractions of inactive 
processes one observes strikingly different behaviors for the reversible and irreversible 
components (Fig. [2j right). The former, in particular, suggests the existence of a second- 
order phase transition at n = 1/(1 + <f>): for n e f < 1 (i.e. in the contracting regime) all 
reversible processes are active, while ipo tre v > strictly in the expanding regime with 
n e f > 1. 

To verify these predictions we have computed optimal growth solutions numerically 
by an extension of the Minover + algorithm defined in (5) that accounts for reversible 
reactions. Details are given in the Appendix. It turns out that the analytic calculation is 
in very good agreement with numerical estimates in the expanding phase (see Fig. [I]). In 
the contracting regime, however, one observes deviations from the analytic curves that 
increase with 0. In the light of the previous discussion, it is reasonable to expect that 
such deviations are related to a breakdown of convexity at g*. To test this hypothesis, 
we have measured, at fixed N and M, the probability II that the linear combination 
of two solutions with coefficients 1/2 is not a solution of the original problem (see Fig. 





Figure 2. Left panel: fraction of intermediate reagents (cq) versus n e f for different 
values of (p. Right panel: fraction of inactive processes (tpo) versus n e / for different 
values of (/). The reversible and irreversible components are shown separately. 




Figure 3. Numerical estimate of the probability II that the linear combination of two 
solutions with coefficients 1/2 doesn't solve the fully reversible Von Neumann problem 
{4> = 1). Values of n — N/M are as reported in the panels. Note that the optimal 
growth rate lies in the contracting (resp. expanding) phase in the left (resp. right) 
panel. System sizes are as follows: N = 10, 25, 50 and 100 increasing in the direction 
of the arrow for n = 1/8 (note that data for N — 10 may suffer from small system 
size); N = 200,500, 1000,2000 and 4000 increasing in the direction of the arrow for 
n = 20. 



|3j). One observes clear signs of convexity breakdown for g < (or p < 1), while the 
solution space appears to be convex both for g > and when g is sufficiently smaller 
than g* . When the optimal solution lies in the contracting regime, II displays strong 
numerical fluctuations close to g* making it hard to identify unambiguously whether a 
single solution survives. Note that, while the transition at g = is correctly located 
by such an analysis, the crossover point seen numerically for g < is just an upper 
bound of the real transition point where the solution space turns from being convex to 
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— theory 
a-B p= 0.985 
s-a p=0.975 
" p=0.95 



3 



4 



S 



Figure 4. Distribution of s for reversible processes for n e f = 1/4 (contracting phase) 
and 4> = 1 (M = 8./V = 800). The analytical prediction (continuous black line) is shown 
together with distributions for increasing values of p < p*. For these simulations, we 
estimate p* ~ 0.986. 



non-convex upon increasing g. 

We have estimated numerically the flux distributions in the contracting phase for 
= 1 and n e j = 1/4 (see Fig. [4]). The non convexity of the solution space makes it 
hard to achieve the optimal growth solutions (especially so for the negative part of the 
flux distribution). In spite of this, our results clearly show that as one gets closer to g* 
the distribution of fluxes develops the double peak structure predicted by the replica 
theory. 

5. Outlook 

Von Neumann's expanding model acquires considerable complexity in the presence 
of reversibility. Our results show in particular that multiple dynamical paths for reaction 
rates exist in contracting phases due to a breakdown of convexity of the solutions' 
space, whereas a convex solution space, with an unique optimal trajectory, characterizes 
expanding phases. Our theory is able to describe the expanding phase and does appear 
to capture some of the salient features of the contracting regime. At the same time, this 
study opens a number of interesting avenues for future research: In first place, it would 
be important to understand whether the non-convex solution space for p < 1 can be 
described by a replica symmetry breaking Ansatz. Secondly, non-convexity poses the 
problem of designing efficient algorithms to sample the solution space in the contracting 
regime. The algorithm presented here is a promising first step in this direction. 
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Appendix. Algorithm to compute p* 



The Minover + algorithm |4] exploits the similarity between Von Neumann's 
conditions and the pattern storage condition in the perceptron to compute p* on any 
graph in the irreversible Von Neumann model where fluxes are semipositive. In this 
case, the procedure to find a solution at a given p is based on iteratively rotating the 
flux vector s in the direction of the least satisfied constraint, similar to |9|: 



(i) At step j = 0: randomly initialize the flux vector s(0); 

(ii) At step j + 1: compute 

/i = arg min c M (j) , cf{j) = V Si(j)(a% - p%) (36) 

(hi) If c Mo (j) > 0, then s(j) is a solution and exit; 
(iv) Else, if c Mo < 0, update fluxes as follows: 

Si (j + 1) = max{0, Si (j) + af o(i) - pb? o{j) } (37) 

and return to (ii). 



If at time j more than one reagent satisfies (36), a single po can be chosen by picking 
one at random with uniform probability among all metabolites having the minimum 
value of & . In addition, at the end of each step it is possible to normalize the flux 
vector appropriately (e.g. Yli s i = N). Iteration of this subroutine for increasing values 
of p allows to compute p* with the desired degree of accuracy. Note that convergence 
to a solution is guaranteed for every p < p* and, in addition, it can be shown that the 
algorithm samples convex solution sets uniformly. 

To extend this procedure to the reversible model it is not sufficient to 



straightforwardly eliminate the lower bound in (37), since the reinforcement term in 
this case is not fixed (as for a standard perceptron) but depends on the flux itself. In 
particular it is given by 

Af = Si e( Si ) « - PK) - s t 9(sm - K) ( 38 ) 
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The non-linearity of the above relation deforms the 'patterns' each time the direction of 
a reaction is reversed, so it is not immediately obvious in which direction it is necessary 
to modify a flux to approach the solution. 

For reversible reactions, we thus introduce a different core algorithm to solve Von 
Neumann's conditions (to be iterated as before over p to calculate p*). Let us set 

Xt( P ) = af-pK , Yt{p)=bt-pat , K{p)=X^p)-Yt{p) (39) 

The first two quantities represent, respectively, the net amount of metabolite p that the 
i-th reaction gives to (if positive) or subtracts from (if negative) the system per unit of 
flux, in straight and reverse direction respectively. The last one, equal to {l+p)(a% — 6f), 
has the same sign of direction that the i-th reaction has to take to produce the largest 
amount of metabolite p per unit of flux. We proceed as follows. 

(i) Step j = 0: randomly initialize the flux vector s(0) (e.g. as a random vector with 
uniformly distributed entries in [—1,1]); 

(ii) Step j + 1: compute 

p = arg min c M (j) (40) 

with 

N 

i=i 

(iii) If > 0, then s(j) is a solution and exit;; 

(iv) Else, if c Mo (j) < 0, update fluxes as follows: 

(a) If min{(Xf°(p),r. W) (p)} < then 

f max{0, Sl (j) + af° O) ~ pK° U) } 
s l {] + 1) = I min{0, Si{j) ~ - P«r 0) )} 

{ Kip) 

(b) Else, if min{(Xf , (p),^°(p)} > 0: 
(b.l) If Si(j)/if°(p) > follow the same instructions as in (a); 
(b.2) Else, if Sj(j)/if°(p) < update fluxes as 

8 i {j + l)=8 i {j) + h?(j>) (43) 
and return to (ii). 

This procedure tries in essence to favor for each reaction the direction allowing for an 
increase of the production of the reagent corresponding to the least satisfied constraint, 
unless both directions are feasible (i.e. unless both Xf° and are positive), in which 
case the flux vector is rotated in the most advantageous direction. Finally, if at the 
previous step the reaction was inactive, the reaction has to be activated in the direction 
that favors maximal production of Pq. 

In a network with partial reversibility, it is possible to update each flux using the 
reversible or the irreversible algorithm according to the type of reaction considered. 
We have been unable to prove convergence (as was instead possible for the Minover + 
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algorithm) and the only support to its effectiveness lies (besides numerous tests) in its 
ability to recover the replica prediction. 



